Unconstrained Variational Determination of the Kohn-Sham Potential 
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Density Functional Theory's Kohn-Sham (KS) potential emerges as the minimizing effective 
potential in an unconstrained variational scheme that does not involve fixing the unknown single- 
electron density. The physical content behind the virtual KS system, that of a non-interacting 
electronic system representing optimally the interacting one, is brought to light. 



I. INTRODUCTION 



Density Functional Theory [1-6] (DFT) has revolu- 
tionised the study of electronic structure. Its apphca- 
tion, especially under the Kohn-Sham (KS) scheme [2,7], 
prevails: in condensed matter physics DFT is, without 
doubt, the method of choice for ab initio calculations; 
in computational chemistry its popularity is broad and 
growing. Still, despite the success and widespread use, 
the physical content behind the virtual KS system is ob- 
scure and its introduction to DFT remains ad hoc. The 
KS system is a system of non-interacting electrons with 
the same ground state (gs) single-particle density as the 
system of interacting electrons of interest. Some ques- 
tions arise: Would one expect, a priori, this virtual sys- 
tem to represent the interacting system faithfully? In 
other words, would one expect, a priori, the KS state to 
be close to the interacting ground state ^E"? (Relative to 
other non- interacting states, in the same way we expect 
the Hartree-Fock (HF) or the Optimised Effective Po- 
tential [8-10] (OFF) state to be close to Further one 
may ask regarding the KS potential: Beyond the well- 
known theorem that the energy of the highest occupied 
orbital is the negative of the ionization energy [11-13], 
would one expect the KS potential to have 'physical re- 
ality' [14] and yield accurate single-particle excitation en- 
ergies [15], as does the OFF potential? [16] 

These questions cause some degree of discomfort, be- 
cause one is inclined to answer negatively (see however 
[15]), yet after several decades of successfully applying 
the KS method, we know that the KS scheme does have 
intuitional value and physical content and it cannot be a 
mere mathematical construct to obtain the single-particle 
density. It is important, and indicative of our level of un- 
derstanding, that there is a way to predict it is so before 
doing the calculation. 

This Letter addresses this issue. It is shown that the 
KS potential is the optimal potential in a minimization 
that does not involve constraining the single-particle den- 
sity. The KS scheme thus arises in a way entirely analo- 
gous to the OEP method. 



II. THEOREM 

To proceed, denote by Ej the j-th eigenstate, 
eigenvalue of the system of N interacting electrons that 
we are interested in (with a non-degenerate gs): 



Ej'fj, J = 1,2,. 



(1) 



where the Hamiltonian which describes the system is 

H^f + f + , (2) 

T is the electronic kinetic energy operator, V°'^ is the 
nuclear Coulomb attraction operator that binds the elec- 
trons and V"^"^ is the electron-electron Coulomb repulsion 
operator. 

Denote also by ^v,j, Ey.j, the j-th eigenstate and 
eigenvalue of a system of N non-interacting electrons 
{y^'^ = 0), confined by a local potential V{t): 



V ^V.j 



where, 



v,j'i>v,j, J 



Hv--T + V. 



1,2, 



(3) 



(4) 



We shall omit the subscript 1 from the gs, '^i — ^, 
^VA ^ ^v, El ^ E, Eva = Ey- The subscript V de- 
notes dependence on V{r) which, can be almost any local 
potential. It is traditional in DFT to restrict the space 
of potentials to be L"^/^ + L°°. Less restricted spaces also 
make sense [4]. 

Reminder of OEP [8,9]: The following inequality holds: 



{^v\H\^v) ~ E > 



(5) 



The optimized effective potential, T/°°p, is that potential 
which minimizes the Ihs of (5) . The non- interacting state 
$yo<!p is expected to be close to ^ (compared with other 
non-interacting states), because ^yoop minimizes 



-E 



(6) 



in a restricted space, while Vl/ is the global minimizing 
state in the whole Hilbert space. Similarly, the Hartree- 
Fock Slater determinant, $hf, which minimizes (6) over 
all Slater determinants, should be close to 5* as well. It 
obviously holds that 



^'OEP > Ehf > E , 



(7) 



1 



where 



EoEP = ($vo,p|iJ[$yoop) , 
EhF = ($HF|i^|$HF) . 



(8) 
(9) 



Another inequaUty also holds, since Ey is the gs energy 



{^\Hv\'if) - Ev > 0. 



(10) 



Ineq. 10 is the dual of (5), as the roles of the Hamilto- 
nians and the states over which we calculate expectation 
values arc reversed. Wc may interpret (10), by thinking 
of "if as being somewhat near the gs of Hy- It makes 
sense to ask about the quality of this approximation and 
to optimize it: Keeping '5 fixed, what is the potential Vg 
which makes the Ihs of (10) as small as possible? 

The gs, , of this optimal potential will again be 
close to ^' and will play a similar role to $yoep or $hf- 
In fact $yoep are related by the duality of the vari- 
ational principles (5, 10) they originate from. They are 
expected to be similarly 'close' to the gs * and the effec- 
tive non-interacting systems they describe to be equiv- 
alent. In a subsequent paper, we demonstrate the ad- 
vantages to employ or $v""»p as a starting point to 
approximate ^' using perturbation theory. 

It is convenient to give a name to the energy difference 
on the Ihs of (10): 



T^[V] = {^\Hv\'^)-Ev 



(11) 



We will assume Vg is not degenerate. The case of a de- 
generate optimal potential Vg will be discussed later on 
and in detail in a subsequent paper. 

Theorem 

The potential Vg which minimizes r*[V] is the Kohn- 
Sham (KS) potential. The ground state of Vg is the KS 
non- interacting state having the same density as ^. 

Proof 

When we vary the potential, V V+e W, its gs changes, 



$1 



<I>y + e $y 



(12) 



The first-order correction $y depends on V^(r) and W{r). 
For the proof, we only need that, for e 
thogonal to $y: 



0, $y IS or- 



(13) 



At the minimum, the first-order variations of T^[T^] in 
(11) must vanish. Upon varying the potential around the 
minimizing value Vg we have: 



lim 



mHvA^)-Ev^]}=0. 



(14) 



Using (13), we obtain 

Therefore, for any VF(r): 

drWir) (p*(r)-py,(r)) =0, 



(15) 



(16) 



where p^{r) is the density of ^' and pv^ (r) the density of 
^Vs- Since the above must hold for any W{r), we have 
for all r: 



p^{r)-pvSr) =0. 



(17) 



The Ihs of the above is the functional derivative of T$ [V] 



at Vg. 



ST^[V] 



6V{r) 



p*(r) -py^(r). 



(18) 



Vs 



which must vanish. 

We have shown that Vg has the same gs density as the 
interacting external potential V™. By definition, the KS 
potential is the non-interacting potential with the same 
gs density as the interacting one [2]. The basic theorem 
of DFT by Hohenberg and Kohn [1] says that we cannot 
have two different potentials (here both non-interacting) 
with the same gs density. Therefore, Vg must be the KS 
potential and the KS state. 



III. DISCUSSION 

The optimal state need not be a single Slater de- 
terminant. Wc may require that it satisfies symmetries 
of H, for example we may choose it an eigenstate of S"^ 
and Sz- 

The value of [V] at the minimum is 



T^[Vg] = {^\T\^/) - ($yjT|$yJ > 0. 



(19) 



One may recognise T>i,[ys] as the kinetic part, Tcl/O*], 
of the correlation energy functional -Ecb*] = (^'|T' + 

The optimal potential Vg is determined with a free- 
dom of a constant. Separately, the energies {'i'lHyJ'^) 
and Ev^ appearing in T5,[ys] do not give any estimate 
of the total energy, E, since they are arbitrary within a 
constant shift. 

The expectation value {^yJH\^Vs) gives an upper 
bound to i^oEP, 



because V°^p minimizes (5). 



(20) 
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In the derivation of (17) we never used that is the gs 
of H. \1/ could be any state and then the minimization of 
T^[V^] would yield a non- interacting potential Vg whose 
gs would have the same single-particle density as 
The question of non-interacting w-representability arises 
here non-trivially: docs the optimal Vg exist in the space 
of potentials considered? Going back to the proof of the 
theorem, rigorously, we do not know whether the min- 
imum of ^^([y] exists for all Only the infimum can 
always be defined and then (17) might not hold exactly. 
However, having an unconstrained variational principle 
(10) to determine the optimal potential, guarantees that 
we may approach numerically the infimum of [V] as 
close as we wish. 

An example will illustrate this: ^! can be chosen an ex- 
cited state of H. Then, the non-interacting potential Vg 
should develop divergences to simulate in its gs density 
the nodes present in the excited state density p^,. We 
may construct a steepest-descent algorithm to minimize 
T^[y]. This algorithm would generate a series of non- 
interacting potentials with spikes of increasing height at 
the nodes of p^i, , so that the gs densities of these poten- 
tials would approach pq, successively closer. 

A similar variational principle to (10) can be used to 
derive the optimal effective spin- up/down potential when 
the interacting system lies in a weak uni-axial magnetic 
field defining the z direction. Here, we still denote the gs 
of the interacting system by Even in the absence of 
an external magnetic field, allowing for different effective 
spin-up/down potentials V {V\ V^) in the energy dif- 
ference T*[F] in (11), in general, will lead after optimiza- 
tion to a lower minimum, thus implying that the non- 
interacting state $(yT yi) offers a better representation 
of 5* than the original state ^v^. Physically, we expect 
this to become relevant when the electrons in the system 
do not form closed shells. The variational principle that 
after minimization will determine the spin-potential is: 

(^IFf^T.yi)!*) - E^y.yt^ > , (21) 

where iJ(yT,vi) is the gs of the non-interacting spin- 
diagonal Hamiltonian 

H(v^yi)=T + V^ + . (22) 

Optimizing separately the Ihs of (21) for the spin-up and 
down potentials, one finds that the optimal spin-potential 
{Vj ,V}) is the KS spin-potential with the correct spin- 
density (unrestricted KS scheme, or KS scheme in Spin- 
DFT [17,18]): 

(p[.,(r),pJ,,(r)) = (p;(r),pi(r)). (23) 

In the absence of an external magnetic field, since the 
minimum of T^[y^,y''^] (Ihs of (21)) is in general lower 
than the minimum of rii,[y], the following inequality 
holds rigorously for the kinetic part of the correlation 
energy functional: 



n[pl,pi]<T,[pl+pi]. (24) 

A similar inequality should hold, at least approximately, 
for the magnitude of the correlation energy functionals. 

We note that a generalised HK theorem holds [19,20] 
for the spin-diagonal Hamiltonian (22), as long as spin- 
up and down potentials are defined with a freedom of a 
spin-constant and the system described by H(^Ytyi-j is 
not perfectly spin-polarized. 

Finally, let us consider the degenerate case. When the 
optimal potential Vs is degenerate, the gs is not 
unique. Then, a non-interacting state in the space of 
degeneracy will have single-particle density equal to p^. 
I do not find this statement very helpful though. Of- 
ten degeneracy is due to open shells and unpaired elec- 
trons. Then, allowing for different spin- up/down effec- 
tive potentials may lift the dc^gcmcracy. If the degen- 
eracy persists even in tlic presence of an effective mag- 
netic field (when Vj 7^ y/), then use of subspacc or 
ensemble KS scheme for degenerate or low excited states 
[21-25] is appropriate. We shall devote another paper 
to the quasi-degenerate case and the development of an 
ab-initio, multi-configurational KS scheme. 

In conclusion, we have shown that the KS system, until 
now an ad hoc construction with rather esoteric and not 
a priori transparent physical content, can be seen as an 
optimal non-interacting system to represent the gs of the 
interacting electronic system, completely analogous or 
dual to the OEP system. This provides a bridge between 
wave-function methods traditionally employed in quan- 
tum chemistry and DFT. For example, using accurate gs 
wave-functions ^ calculated by other methods and mini- 
mizing the energy difference functional [V] , we can ob- 
tain accurate KS potentials for the studied systems and 
assess the accuracy of existing approximate exchange- 
correlation potentials [26 29]. More importantly, the 
prospect to explore and devise new expressions for the 
correlation energy and potential based on the variational 
principle (10) appears promising. 

K. Burke, R. van Leeuwen, M.W. Long and C.A. Ull- 
rich communicated helpful comments on an early version 
of this manuscript. 
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